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^ . The basics of correlation femtoscopy, recent results from femtoscopy in relativistic heavy 

O ! ion collisions and their consequences are shortly reviewed. 
(N 
O 

O ' 1- Introduction 
^ I 

^ ■ The momentum correlations of two or more particles at small relative momenta in their 
center-of-mass (cm.) system are widely used to study space-time characteristics of the 
production processes on a level of fm = 10~^^ m, so serving as a correlation femtoscopy 



tool. Particularly, for non-interacting identical particles, like photons or, to some extent, 
^ I pions, these correlations result from the interference of the production amplitudes due to 
^ ' the symmetrization requirement of quantum statistics (QS) |H2|3|4|5] . 

There exists P^3,.4...5,.6-7i|l8j an analogy of the momentum QS correlations of photons 
^ . with the space-time correlations of the intensities of classical electromagnetic fields used 
in astronomy to measure the angular radii of stellar objects based on the superposition 
principle (HBT effect) j^]. This analogy is sometimes misunderstood and the momentum 
correlations are mixed up with the space-time (HBT) correlations in spite that their 
orthogonal character and thus the absence of the former in astronomy measurements 
due to extremely large space-time extent of stellar objects (and vice versa) was already 
pointed out in early paper [H]. Note that though the space-time (HBT) correlations 
are absent in the subatomic measurements, they can still be used in laboratory as an 
intensity-correlation spectroscopy tool (see jl(J|llj and references therein). 

The momentum QS correlations were first observed as an enhanced production of the 
pairs of identical pions with small opening angles (GGLP effect |1]). Later on, Kopylov 
and Podgoretsky settled the basics of correlation femtoscopy in more than 20 papers (see 
a review [Hj) and developed it as a practical tool; particularly, they suggested to study the 
interference effect in terms of the correlation function, proposed the mixing techniques 
to construct the uncorrelated reference sample and clarified the role of the space-time 
characteristics of particle production in various physical situations. 
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The momentum correlations of particles emitted at nuclear distances are also influ- 
enced by the effect of final state interaction (FSI) |12|13|14|15|16j . Thus the effect of 
the Coulomb interaction dominates the correlations of charged particles at very small 
relative momenta (of the order of the inverse Bohr radius \a\^^ of the two-particle sys- 
tem), respectively suppressing or enhancing the production of particles with like or unlike 
charges. Though the FSI effect complicates the correlation analysis, it is an important 
source of information allowing for the coalescence femtoscopy (see, e.g., |17ll8ll9'PDj ). the 
correlation femtoscopy with unlike particles |14|15llb'] including the access to the relative 
space-time asymmetries in particle production |21l22l2;-il24l25l2eil27l2SI29i;-i()i;-ili;-i2j and 
a study of strong interactions between specific particles |28|31|32j . 

In this review, I will concentrate on the assumptions behind the correlation femtoscopy 
formalism and discuss the recent results obtained from the femtoscopy analysis of like 
and unlike particle correlations in relativistic heavy ion collisions. One can inspect recent 
reviews |32|33|34j for a number of important topics that are not touched here, such as 
non-Gaussian tails, imaging techniques, correlations of penetrating probes, comparison of 
different colliding systems or spin correlations. 

2. Assumptions 

The two-particle correlation function TZ{pi,p2) is usually defined as a ratio of the mea- 
sured two-particle distribution to the reference one obtained by mixing particles from 
different events of a given class, normalized to unity at sufficiently large relative mo- 
menta. The space-time information contained in the momentum correlations is usually 
extracted based on the following assumptions: 

(i) The mean freeze-out phase space density (/) is assumed sufficiently small so that 
only the mutual QS and FSI effects can be considered when calculating the correlation 
function of two particles emitted with a small relative momentum Q = 2k* in their cm. 
system. This two-particle approximation may not be justified for the rare pairs associated 
with large phase-space density fluctuations and also in low energy heavy ion reactions 
when the particles are produced in a mean field of residual nuclei. To deal with this field 
a quantum adiabatic (factorisation) approach j35|36j or the transport simulations (see, 
e.g., |37|38|39j ) can be used. At high energies, there are also some attempts to account 
for the mean field effects in eikonal |4U|41|42] and optical potential |43|44j approaches, 
the latter giving some tens of percent change of the outgoing phase space density profiles 
in heavy ion collisions at RHIC assuming however unrealistic in-medium masses at the 
kinetic freeze-out and neglecting a rapid time dependence of the potential. As for (/), it 
increases with energy and for central lead-lead or gold-gold collisions seems to saturate 
at the highest SPS energy |33|45j (see, however, [lEI)- The saturated (/) is substantially 
smaller than unity for pions with pt > 0.2 Gev/c so pointing to negligible multiboson 
effects (see, e.g., |47|48j ) in this pt-iegion. 

(ii) The momentum dependence of the one-particle emission probabilities is assumed 
inessential when varying the particle four-momenta pi and p2 by the amount character- 
istic for the correlations due to QS and FSI. This smoothness assumption, requiring the 
components of the mean space-time distance between particle emitters much larger than 
those of the space-time extent of the emitters, is well justified for heavy ion collisions. 
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(iii) An independent or incoherent particle emission is assumed. This assumption is 
quite reasonable for a dominant part of particle pairs produced in heavy ion collisions and 
is consistent with the observed strength of two- and three-pion correlation functions (see, 



(iv) To simplify the calculation of the FSI effect, the Bethe-Salpeter amphtude de- 
scribing two particles emitted at space-time points Xj = {ti,rj} and detected with four- 
momenta Pi is usually calculated at equal emission times in the pair cm. system; i.e. 
the reduced non-symmetrized Bethe-Salpeter amplitude (with the removed unimportant 
phase factor due to the cm. motion), depending only on the relative four-coordinate 
Ax ~ Xi — X2 = {t,r} and the generalized relative momentum q = q — P{qP)/P^ 
{q = Pi — P2, P = Pi + P2 and qP = mi^ — m2^; in the two-particle cm. system, P = 0, 
q = {0,2k*} and Ax = {t*,r*}), is substituted by a stationary solution 'il:_\^J{r*) of the 
scattering problem having at large distances r* the asymptotic form of a superposition of 
the plane and outgoing spherical waves (the minus sign of the vector k* corresponds to 
the reverse in time direction of the emission process). This equal time approximation is 
valid on conditions |14|52j \t*\ -C m2,ir*^ for sign(t*) = ±1 respectively. These conditions 
are usually satisfied for heavy particles like kaons or nucleons. But even for pions, the 
t* = approximation merely leads to a slight overestimation (typically < 5%) of the 
strong FSI effect [H2] and, it doesn't influence the leading zero-distance (r* -C |a|) effect 
of the Coulomb FSI. 

Then, for non-identical particles. 



s 

similar to the Fermi factor taking into account the Coulomb FSI in /3— decay. For identical 
particles, the amplitude in Eq. (Q) enters in a symmetrized form: 



The averaging in Eq. is done over the four-coordinates of the emitters at a given total 
spin S of the two-particles, ps is the corresponding population probability, J2sPs = 1- For 
unpolarized particles with spins si and S2 the probability ps = (25'+l)/[(2si + l)(2s2 + l)]. 
Note that in the case of small k*, we are interested in, the short-range interaction is 
dominated by central forces and s- waves so that, neglecting a weak spin dependence of 
the Coulomb interaction, the spin dependence of the two-particle amplitude enters only 
through the total spin S. 

(v) The two-particle approximation in (i) and FSI separation in the Bethe-Salpeter 
amplitudes of the elastic transitions 1-1-2^1 + 2 implies a long FSI time as compared 
with the characteristic production time, i.e. the channel momentum k* much less than 
typical production momentum transfer of hundreds MeV/c. In fact, the long-time FSI 
can be separated also in the inelastic transitions, 1 + 2 3 + 4, characterized by a 
slow relative motion in both entrance and exit channels |53|52j . The necessary condition 
is an approximate equality of the sums of particle masses in the channels a = 1 + 2 
and /3 = 3 + 4. In the presence of such transitions the amplitudes in Eqs. (0) and Q 
should be substituted by the solutions of the two-channel scattering problem, i.e., formally, 
S S,aa and S,l3a. In practice, the particles 1,3 and 2,4 are members of the same 



e.g., mmm). 



^(pi,P2) = EP5(|Clt\r*)r)5, 



(1) 



(2) 
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isomultiplets (as, e.g., in the transition n p ^ vr^n) so that one can put ps,aa = Ps,i3a 
and assume the same r*-distributions in the channels a and /3. 

3. Femtoscopy with identical particles 

For identical pions or kaons, the effect of the strong FSI is usually small and the effect of 
the Coulomb FSI can be in first approximation simply corrected for (see [SI] and references 
therein). The corrected correlation function is determined by the QS symmetrization only 
(see Eq. (j2I) and substitute the non-symmetrized amplitude by the plane wave e*''^/^): 



Its characteristic feature is the presence of the interference maximum at small components 
of the relative four-momentum q with the width reflecting the inverse space-time extent 
of the effective production region. 

Note that the on-shell constraint go = vq (v = P/Pq is the pair velocity) makes the 
g-dependence of the correlation function essentially three-dimensional ^ (particularly, 
in pair cm. system, qAx = — 2k*r*) and thus makes impossible the unique Fourier 
reconstruction of the space-time characteristics of the emission process. However, within 
realistic models, the directional and velocity dependence of the correlation function can be 
used to determine both the duration of the emission and the form of the emission region 
|2|3|4|5j . as well as - to reveal the details of the production dynamics (such as collective 
flows; see, e.g., |55l56l57j and reviews |58|59j ). For this, the correlation functions can be 
analyzed in terms of the out (x), side (y) and longitudinal (z) components of the relative 
momentum vector q = {qx,%,(lz} |6m61|62|63j : the out and side denote the transverse 
components of the vector q, the out direction is parallel to the transverse component of the 
pair three-momentum. The corresponding correlation widths are usually parameterized 
in terms of the Gaussian correlation radii Ri, 



and their dependence on pair rapidity and transverse momentum is studied. The corre- 
lation strength parameter A represents a fraction of the pairs of identical pions or kaons 
emitted by independent short-lived sources. The form of Eq. (j^ assumes azimuthal sym- 
metry of the production process |58|(j(Jj . Generally, e.g., in case of the correlation analysis 
with respect to the reaction plane, all three cross terms q^qj contribute |64j . 

It is well known that particle correlations at high energies usually measure only a small 
part of the space-time emission volume, being only slightly sensitive to its increase related 
to the fast longitudinal motion of particle sources. In fact, due to limited source decay 
momenta of few hundred MeV/c, the correlated particles with nearby velocities are emit- 
ted by almost comoving sources and so - at nearby space-time points. The dynamical 
examples are sources-resonances, colour strings or hydrodynamic expansion. To substan- 
tially eliminate the effect of the longitudinal motion, the correlations can be analyzed in 
terms of the invariant variable Q = {—(fY^'^ = 2k* and the components of the momentum 
difference in pair cm. system (q* = Q = 2k*) or in the longitudinally comoving system 
(LCMS) jnS]- In LCMS each pair is emitted transverse to the reaction axis so that the 
generalized relative momentum q coincides with q* except for the component q^ = jtcix^ 
where 7* is the LCMS Lorentz factor of the pair. 



T^{Pi,P2) = 1 + (cos(gAx)). 



(3) 



^(Pi,P2) = 1 + Xexp{-Rlql - Rlql - R^q^ - 2Rl^qxqz) 



(4) 
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Particularly, in the case of one-dimensional boost invariant expansion, the longitudinal 
correlation radius in the LCMS reads [57j Rz ~ {T /mtY^'^T, where T is the freeze-out 
temperature, r is the proper freeze-out time and rrit is the transverse particle mass. In 
this model, the side radius measures the transverse radius of the system while the square 
of the out radius gets an additional contribution (pf/mt)^Ar^ due to the finite emission 
duration At. The additional transverse expansion leads to a slight modification of the 
Pt-dependence of the longitudinal radius and - to a noticeable decrease of the side radius 
and the spatial part of the out radius with pt. Since the freeze-out temperature and the 
transverse flow determine also the shapes of the mt-spectra, the simultaneous analysis of 
correlations and single particle spectra for various particle species allows to disentangle 
all the freeze-out characteristics (see a review [5S]). It appears that with the increasing 
energy of heavy ion collisions from AGS and SPS up to the highest energies at RHIC, 
the correlation data show rather weak energy dependence |66|67j and point to the kinetic 
freeze-out temperature somewhat below the pion mass, a strong transverse flow (with 
the mean transverse flow velocity at RHIC exceeding half the velocity of light), a short 
evolution time of 8-10 fm/c and a very short emission duration of about 2-3 fm/c (see, 
e.g., a recent review [SH]). The short evolution and emission duration at RHIC are also 
supported by the correlation analysis with respect to the reaction plane |.67 68 . 

The small time scales at RHIC were not expected in traditional transport and hy- 
drodynamic models (see, e.g., |69|7Uj ) and may indicate an explosive character of particle 
production |7H72j or, large parton cross sections and early hadronization jTlj as in the 
corresponding successful multi-phase and hadron transport models. In view of a successful 
hydrodynamical description of the elliptic flow at RHIC (see, however, |j75j) the failure of 
hydrodynamics to explain the femtoscopy data is often considered as a puzzle. However, 
it is still possible that the overestimation of the outward and longitudinal radii will be 
cured in the full three-dimensional hydrodynamic calculations with the modified initial 
conditions (earlier thermalization time or initial transverse flow |76|77j ). tuned equation 
of state with the account for chemical non-equilibrium after hadronization |7til77t78,.79j 
and a better modeling of the freeze-out process (including continuous emission during the 
hydrodynamic stage |8U|81j and coupling to the hadron cascade |82|83j ). This hope is 
supported by a better description of the correlation radii in a recent three-dimensional 
hydrodynamic model |M1 and by a successful description of the SPS and RHIC data in 
a number of hydro- motivated parametrizations (see a brief review jHH])- The success of 
these simple parametrizations can be caused by the formation of particle spectra and 
correlation radii during an early stage of the kinetic freeze-out (see, e.g., |83|73|74j ). It 
was argued [SI] that such an early spectra formation can be related with the fact that the 
solution of the non-relativistic Boltzmann equation at spherically symmetric initial con- 
ditions exactly coincides with free streaming (see, e.g., [86j). Indeed, numerical solutions 
of relativistic Boltzmann equation at anisotropic initial conditions showed only a slight 
deviation from free streaming even in the case of a large number of rescatterings [HZ] . 

4. Femtoscopy with unlike particles 

The complicated dynamics of particle production, including resonance decays and parti- 
cle rescatterings, leads to essentially non-Gaussian tail of the r*-distribution. Therefore, 
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due to different r*-sensitivity of the QS, strong and Coulomb FSI effects, one has to be 
careful when analyzing the correlation functions in terms of simple models. Thus, the 
QS and strong FSI effects are influenced by the r*-tail mainly through the correlation 
strength parameter A while, the shape of the Coulomb FSI is sensitive to the distances 
as large as the pair Bohr radius \a\ (hundreds fm for the pairs containing pions). These 
problems can be at least partially overcome with the help of transport code simulations 
accounting for the dynamical evolution of the emission process and providing the phase 
space information required to calculate the QS and FSI effects on the correlation function. 

Thus, in a preliminary analysis of the NA49 correlation data from central Pb + Pb 
158 AGeV collisions 28^, the transport RQMD v.2.3 code was used. To account for a 
possible mismatch in (r*), the correlation functions were calculated with the space-time 
coordinates of the emission points scaled by 0.7, 0.8 and 1. The scale parameter was then 
fitted using the quadratic interpolation. The fits of the vr+vr^, ir^p and ir'p correlation 
function indicate that RQMD overestimates the distances r* by 10-20% thus indicating 
an underestimation of the collective flow in this model. 

Recently, there appeared data on pA correlation functions from Au + Au experiments 
E985 at AGS [Hi, STAR at RHIC ^ and Pb + Pb experiment NA49 at SPS CERN 
PU] . As the Coulomb FSI is absent in pA system, one avoids here the problem of its 
sensitivity to the r*-tail. Also, the absence of the Coulomb suppression of small relative 
momenta makes this system more sensitive to the radius parameters as compared with 
pp correlations [Sni- In spite of rather large statistical errors, a significant enhancement is 
seen at low relative momentum, consistent with the known singlet and triplet pA s-wave 
scattering lengths. In fact, the fits using the analytical expression for the correlation 
function ^3] yield the correlation radii of 3-4 fm in agreement with the radii obtained 
from pp correlations in the same experiments. These radii are smaller than those obtained 
from two-pion and two-kaon correlation functions at the same transverse momenta and 
are in qualitative agreement with the approximate rrit scaling expected in the case of the 
collective expansion. 

5. Accessing strong interaction 

In case of a poor knowledge of the strong interaction, which is the case for meson- 
meson, meson-hyperon or hyperon-hyperon systems, it can be improved with the help of 
correlation measurements. 

In heavy ion collisions, the effective radius tq of the emission region can be considered 
much larger than the range of the strong interaction potential. The FSI contribution to 
the correlation function is then independent of the actual potential form |14|90j . At small 
Q = 2k*, it is determined by the s-wave scattering amplitudes f^{k*) at a given total spin 
5* ^Ij. In case of \f^\ > tq, this contribution is of the order of |/'^/to|^ and dominates 
over the effect of QS. In the opposite case, the sensitivity of the correlation function to 
the scattering amplitude is determined by the linear term /tq. 

The possibility of the correlation measurement of the scattering amplitudes has been 
demonstrated |^ in a preliminary analysis of the NA49 7r"'"7r~ correlation data within 
the RQMD model. The fitted strong interaction scale, redefining the original s-wave 
scattering length of 0.23 fm, appeared to be significantly lower than unity: 0.63 ± 0.08. 
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To a similar shift (~ 20%) point also the recent BNL data on K14 decays [HI]. These 
results are in agreement with the two-loop calculation in the chiral perturbation theory 
with a standard value of the quark condensate [22] • 

Recently, also the singlet AA s-wave scattering length has been estimated |28l29j based 
on the fits of the NA49 data from Pb + Pb collisions at 158 AGeV. Though the fit results 
are not very restrictive, they likely exclude the possibility of a large singlet scattering 
length comparable to that of ~20 fm for the two-nucleon system. 

The STAR experiment measured for the first time the pA and pA correlation functions 
and performed simultaneous fit of the correlation radius and the spin-averaged s-wave 
scattering length. The fitted imaginary part of the scattering length of ~ 1 fm is in 
agreement with the pp results (thus pointing to about the same pA and pp annihilation 
cross sections) while the real part appears to be more negative ^Tj . 

6. Accessing relative space-time asymmetries 

The correlation function of two non-identical particles, compared with the identical 
ones, contains a principally new piece of information on the relative space-time asymme- 
tries in particle emission j21j . Since this information enters in the two-particle amplitude 
through the terms odd in k*r* = pl{rl — r^), it can be accessed studying the 
correlation functions TZ+i and TZ-i with positive and negative projection k* on a given 
direction i or, - the ratio IZ^i/lZ^i. For example, i can be the direction of the pair veloc- 
ity or, any of the out (x), side (y), longitudinal (z) directions. In LCMS, r* = rj except 
for r* = Ax* = 7t(Aa; — vtAt), where jt and Vt are the pair LCMS Lorentz factor and 
velocity. One may see that the asymmetry in the out (x) direction depends on both space 
and time asymmetries (Ax) and (At). In case of a dominant Coulomb FSI, the intercept 
of the correlation function ratio is directly related with the asymmetry (r*) scaled by the 
Bohr radius a = {^ziZ2e^Y^: TZj^i/TZ^i ^ 1 + 2{r*)/a. 

A review of the simulation studies of the method sensitivity and the experimental results 
can be found elsewhere |32j. Here we only note that the out correlation asymmetries 
between pions, kaons and protons observed in heavy ion collisions at CERN SPS and 
BNL RHIC are in agreement with practically charge independent meson production and 
a negative (Ax) and/or positive c{At) on the level of several fm (assuming mi < 7712) 
|28|29|3(J|3H32j . In fact they are in quantitative agreement with the RQMD transport 
model as well as with the hydro-motivated blast wave parametrization, both predicting 
the dominance of the spatial part of the asymmetries generated by large transverse flows. 
Particularly, the RQMD simulation of central Pb + Pb collisions in conditions of the 
experiment NA49 at 158 AGeV yields practically zero asymmetries for tt+tt^ system 
while, for vr^p systems, (Ax) = —5.2 fm, (At) = 2.9 fm/c, (Ax*) = — 8.5/m. Besides, it 
predicts (x) increasing with particle pt or Vt = pt/rrit, starting from zero due to kinematic 
reasons. The asymmetry arises because of a faster increase with vt for heavier particle. In 
fact, the hierarchy (x^r) < (xk) < (xp) is a signal of a universal transversal collective flow 
[2H|; one should simply take into account that the mean thermal velocity is smaller for 
heavier particle and thus washes out the positive shift due to the flow to a lesser extent. 

These conclusions can be expressed in a simple analytical form using a hydro-motivated 
parametrization. Thus, assuming the longitudinal-boost invariance, a linear non-relativistic 
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transversal flow velocity profile jSp = Poft/f'o, the local thermal momentum distribution 
characterized by the kinetic freeze-out temperature T and the Gaussian density profile 
exp{—r^/{2rl)), one confirms a faster rise of (x) with Vf for heavier particles j|32|931: 



The size of the shift (x) and its mass dependence are weakened when introducing the 
transverse temperature gradient (AT/T)^ like in the Buda-Lund model (this gradient 
just adds to the denominator in Eq. (0)) [H^- The observation of large correlation 
asymmetries at RHIC may thus point against the existence of a hot emitting central zone 
and the corresponding gradient (AT/T)^ ~ 1 indicated by the analysis of RHIC data 
within the Buda-Lund model [HHj. Further quantitative studies are needed to clarify this 
question. 

7. Conclusions 

Wealth of data on momentum correlations of various particle species (vr^, K^^, p^, A, S) 
is available and gives unique space-time information on production characteristics includ- 
ing collective flows. Rather direct evidence for a strong transverse flow in heavy ion 
collisions at SPS and RHIC is coming from unlike particle correlation asymmetries. Be- 
ing sensitive to relative time delays and collective flows, the correlation asymmetries can 
be especially useful to study the effects of phase transitions. Weak energy dependence of 
correlation radii contradicts to traditional hydrodynamic calculations which overestimate 
longitudinal and, especially, outward radii at RHIC. It remains to be clarified whether 
this can be cured in full three-dimensional hydrodynamic calculations with a more refined 
treatment of initial conditions, equation of state and freeze-out process. A number of 
succesful hydro motivated parametrizations give useful hints in this direction. The mo- 
mentum correlations between specific particles yield also a valuable information on the 
strong interaction hardly accessible by other means. 
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